WEYL ASYMPTOTICS FOR TENSOR PRODUCTS OF OPERATORS 

AND DIRICHLET DIVISORS 
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Abstract. We study the counting function of the eigenvalues for tensor products of oper- 
ators, and their perturbations, in the context of Shubin classes and closed manifolds. We 
emphasize connections with problems of analytic number theory, concerning in particular 
generalized Dirichlet divisor functions. 



1. Introduction 

As well known, there are deep connections between spectral theory and analytic number 
theory. One main topic is given by Weyl formula for self-adjoint partial differential operators 
or pseudo-differential operators. Namely, the leading term in the expansion of the counting 
function N(X) of the the eigenvalues < A is recognized to be proportional to the volume of 
the region defined by the A-level surfaces of the symbol, and in turn to the number of the 
lattice points belonging to the region. Even, for relevant classes of operators, each point of 
the lattice corresponds exactly to one of eigenvalues, counted according to the multiplicity, 
and the computation of N(X) leads in a natural way to problems of number theory. Let us 
refer for example to JSJE], |12j-|22|. [28 -[32j. In this order of ideas, the attention will be 
fixed here on operators of the form of tensor products 

(1.1) P = P 1 ®...®P P 

where the operators Pj, j = 1, . . . ,p, are self-adjoint, say strictly positive (pseudo-differential) 
operators on corresponding Hilbert spaces with eigenvalues {A^}^ =1 ,j = I,..., p. Then, the 
eigenvalues of P are products of the form A^.-.A^ and the eigenfunctions are tensor products 
of the corresponding eigenfunctions, cf. [8j |29] for the general functional analytic setting. 
Hence, 

(1.2) N P (\) = #{(k ll ...,k p )GW: \^\%...\^ <\}. 

The computation of N p (X) meets then some classical divisor counting problems. To give a 
simple example, consider the Hermite operators 

(1-3) H j = ±(-% j +xj) + ~, 3 = 1,2. 

Writing for short H\ and H2 for Hi <S> I2 and I\ ® H2, we define the tensorized Hermite 
operator H = H\ <g> H2. In applications, H is sometimes used as a substitute for the standard 
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two dimensional Hermite operator Hi + H2, producing the same eigenfunctions, i.e., two 
dimensional Hermite functions. The distribution of eigenfunctions, counted with multiplicity, 
is however quite different, being related to the distribution of the prime numbers. In fact, the 
eigenvalues of the one-dimensional Hermite operator, normalized as above, are the positive 
integers; therefore, (|l,2p reads in this case as 

[A] 

(1.4) N H (X) = D(X) = Y,d(n), A > 1, 

71=1 

where d(n) denotes the number of divisors of n and [A] stands for the integral part of A. 
Dirichlet proved in 1849 that 

(1.5) D(X) = AlogA + (2 7 -l)A + £(A), 

where 7 is the Euler-Mascheroni constant and E(X) = 0{X 1 / 2 ). The first term on the right 
hand side of fjl .5j) can be easily recognized as the volume of the hyperbolic region defined by 
the symbol of H, whereas the optimal growth order of the rest E(X) is a long-standing open 
problem in the analytic theory of numbers, see for example [Tj [161 EQl EH 133] . 

Natural generalizations of the Hermite operators Hj in fjl .3|) are the global pseudo-differen- 
tial operators of M. Shubin [171 El [251 [32] . If Pj 1S globally elliptic self-adjoint in these classes, 
then the Weyl formula yields 

(1.6) N Pj (X) ~ AjX a i , 

where ctj = 2nj/rrij, with mj the order of Pj and rij the space dimension. The constant 
Aj depends on the symbol of Pj, according to the Weyl formula. Note that the tensorized 
product in (jl.ip is not any longer globally elliptic on W 1 , n = n\ + ... + n p . 

The first aim of the present paper will be to deduce from (|1 .6j) an asymptotic expansion for 
the spectral counting function Np(X) in ()1.2j) . Particular attention will be devoted to lower 
order terms of the asymptotic expansion for some particular cases. As an example, define Hj 



as in (jl.3p and consider now 

(1.7) H^ = H^H^, 

where (3 = , ^2) is a couple of positive integers with P2 ^ Pi- Then, we shall prove that 

(1.8) N H0 (X) = C(/V/3i)A 1//31 + C(/3i//3 2 )A 1 / & + 0(X^ + ^), 



where Q{z) is the Riemann zeta function, analytically continued in the complex plane for 
z + \. 

As we shall also detail in the paper, parallel results can be obtained when Pj in (jl.ip are 
elliptic self-adjoint pseudo-differential operators on a smooth compact manifold. In this case 
(II. 6p is valid with aj = rij/rrij, see [15] , 

Let us finally describe what, to the best of our knowledge, was already known about tensor 
products of pseudo-differential operators and their spectrum, as well as what is new in our 
paper. An algebra of "bisingular" pseudo-differential operators on the product of two mani- 
folds Mi x M2, containing Pi®P2 with Pi and P2 beign classical pseudo-differential operators 
on Mi,M2, respectively, was studied by Rodino [28J in connection with the multiplicative 
property of the Atiyah-Singer index pQ. The spectral properties of this class were recently 
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studied by Battisti [5j. The variant for the Shubin type operators has been considered in [6]. 
These results give a general framework for the study of the example (|1.4p with the expansion 
(jl.5p . and provide as well the leading term in the expansion (| 1 . 8 1) for the example (jl.7p . Let 
us also mention the articles |144 I15j . where starting from the twisted Laplacian of M. W. 
Wong [34J, similar problems of Dirichlet divisor-type were met. The operators in [T5J[T5l[3l] 
are not tensor products, but they can be reduced to the form (jl.ip by conjugating with a 
Fourier integral operator, cf. [13]. 

From the point of view of Mathematical Physics, Kaplitskii [22] has independently studied 
the spectral properties of operators on the torus T 2 with principal part 



obtaining for the counting function estimates of type (jl.5p . Reference therein is made to 
Arnold [3], suggesting to transfer the Weyl formula to hyperbolic equations. The results in 
|22j can be essentially regarded as a particular case of those from [5]. Expansions of the type 
(jl.5p appear also in the recent paper of Coriasco and Maniccia [TO] concerning the spectrum 
of the so-called SG-operators. 

Summing up, the results mentioned above cover the case of products of two operators, 
P = Pi ® P2, except for the computation of lower order terms in the expansions, cf. (jl.8p . 
Thus, our attention will be mainly focus on the case p > 3 of (jl.ip and lower order terms. 

In the present paper, the attention will be rather addressed to results of (elementary) 
analytic number theory, which we shall present in Section [2] in detail; they are new by them- 
selves, we believe. The applications to spectral theory will be given in the conclusive Section 
[3j We shall not construct here an algebra of (multisingular) pseudo-differential operators 
containing Pi (g) ... <g) P p for p > 3. Computations are cumbersome, involving a stratified 
calculus of the type of that from |24[ [2B| [30], 131] . occurring in other contexts. Finally, we 
shall limit ourselves instead to consider perturbations of the type P + Q, where Q is a lower 
order pseudo-differential operator. 



We study in this preparatory section the asymptotic behavior of some counting functions 
of "multi-divisor" type. They will be very helpful when applied to spectral asymptotics of 
various examples of "multi-singular" operators. 

2.1. Counting functions of products of sequences. We start by considering the follow- 
ing general question. Let {X^}^ =1 , j = l,...,p, be non-decreasing sequences of positive 
real numbers. The sequences are rather arbitrary and they are not necessarily linked to any 
operator. 

Assuming that we have some knowledge about each of the counting functions 



P = P x ®Py=dl 



2. Asymptotics of some counting functions 



(2.1) 
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we would like to obtain asymptotic information about the counting function of the p products 
of the elements of the sequences, namely, 

(2.2) N(X):= Y, 1 = #{(*!,...,**) 6 N*: A^A^ . . . < a} . 

a«a1 2) ...a^<a 

The next simple proposition tells us that it is always possible to find the asymptotic 
behavior of (|2.2p whenever there is a block of counting functions (|2.ip with dominating 
asymptotic behavior. 

Proposition 2.1. Suppose that there are non-negative numbers r < a and indices j\, . . . ,j u , 
where 1 < v < p, such that 

(2.3) N jq (X)~ A jq X a , A^oc, g=l,. ..,«/, 
with Aj q 7^ 0, and 

(2.4) Nj(X) = 0(X T ), A^oo, j£{ji,...,ju}. 
Then, the counting function (I2.2p has asymptotic behavior 



(2.5) N(X) ~ AX a { ^ , A -> oo, 



We will divide the proof of Proposition 12.11 into two lemmas. The first lemma deals with 
the case in which all counting functions have asymptotic behavior of the same order. 

Lemma 2.2. If Nj(X) ~ AjX a , with a > and Aj ^ 0, for j = 1,2, . . . ,p, then §Z3§ has 
asymptotics 



N(X)~AX a { S % , A^oo, 
(p- 1)! 



where A = nf=i A? • 



Proof. We proceed by induction. Assume that 

w ^ (p-2 ! 

a^a< 2 )...a(- 1 »<a ^ 7 
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with A = Yl^Zi Aj. We then have, 

k (p) 

A<*><A 

AaP- 2 



N(X) = £ N(X/X^) 



p ^T E (A/A?V(log(A/A^)r 2 + £ o((A/A?Vlog p - 2 A) 



y — 



+ 0(A-Iogr a A)- E ^ 



(p-2)! 



^ (A/Af ) a (log(A/A^)r 2 + o(A a log* 3 - 1 A) + 0(A Q log^ 2 A) 



A-<A 



; „ - 2)7 I (V*) a (log(A/t)) p - 2 aW p (t) + o(A Q log*- 1 A) 
A" j\log(\/t)r 2 ^>dt + (A« log?' 1 A) 

g ( P T 2 ) (-l)^(log A)- - [ **$L* + o(A« log-' A) 
(alogATi ^ - 2\ (-If 

(P-2)! 7o (P-1)! 



We also have, 
Lemma 2.3. If 

Mi(A)= Yj l = 0(A r ) and M 2 (A) = l~£A a log 6 A, A -> oo 

where < t < a, B ^ 0, and b > 0, then 

M(A) = ^ l~SA a log 6 A, A^oo, 

Proof. Observe first that 

y^L_ = r t - a dM l (t)=a r^-dt 
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is convergent because t~ l ~ a M\(t) = 0(t~ 1 ~^ a ~ T ^). Now, 
M(A) 



£ M 2 (X/^) = BX a log b X-B H l ^± d M l {t) + o{\ a \og b \) 

= BX a log 6 A + O (V(log 6 A) £ -^^j + o(X a log 6 A) 
- BX a log b A, 

as claimed. □ 

We can now prove Proposition 12.11 

Proof of Proposition \2. 1\ Let • • • , j p } = {1,2, . . . ,p} \ {ji, . . . ,j u } We arrange the two 

sequences of products 

m n ^ f[x^ 

q=u+l q=l 

in two non-decreasing sequences {/A an d {^^IfcLi; respectively, where each element in 

these sequences is repeated as many times as it can be represented as in (|2.7p . The hypothesis 
([23D and Lemma E2] yield 

M 2 (A) = l~B\ a \og'- 1 \, A^oo, 

where B = (a u ~ 1 j(y — 1)!) Y\q=i A?V On ^ ne other hand, using (12. 4p . one easily verifies that 

Mi (A) = 1 = 0{X T log p - y A), A ->■ oo. 

Applying Lemma 12.31 and noticing that 

£*=( n (e— - 

we obtain the asymptotic formula ()2.5f) with the constant (j2.6|) . □ 

2.2. Remainders. We now study the remainder in (I2.5p . We impose stronger assumptions 
than (12. 3D on the leading counting functions. 

We start by looking at the case when a single counting function dominates all others. 

Proposition 2.4. Assume that there are non-negative numbers t < i] < a and an index 
I £ {1, ... ,p} such that Nj(X) = 0(A T ), for j ^ /, and Ni satisfies 

(2.8) N l (X) = A l X a + 0(X^), A^oo, 
with A\ ^ 0. Then, 

(2.9) iV(A) = ylA Q + 0(X V ), A — ^ oo, 
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where 



/ 00 1 

(2.10) A = A t H i^j- 



&\ \ k =i 

Proof. By renaming the sequences, we may assume that 1 = 1. We use a recursive argument. 
Suppose that we have already established 

N(X) = l = AX a + 0(X^ 

aW'..^- 1 ' 

fc l fc 2 k p— 1 

Since for any b > r 



we have 



^(a?) 6 Ab A tb+1 

N(X)= N{X/X (p) ^ 



k 



A«<A 



\&=i l^fc J / \(p)<-\ 



= ^4A ai +0(X T ) + 0(X 71 ), 
which shows (I2~9i □ 



For the analysis of the remaining case, we will employ a complex Tauberian theorem of 
Aramaki [2]. 

Proposition 2.5. Assume there are non-negative numbers t < a and indices ji,...,j u , 
where 1 < v < p, such that 

(2.11) N ]q (X)=A Jq X a + 0(X T ), A — )► 00, ? = l,...,i/, 

with Aj q 7^ 0, and p?.^| ) holds. Then, there exists rj < a such that (I2.2j) has asymptotic 
expansion 

+ 0(A r < 

! (I — I I 1 \ (I ~ I \ ' I 

(alogA)^ 1 ^ 2 



(2 ' 12) ^ = S(^5f(s)"(t) 



g=0 

where A is given by (|2.6p and i/ie constants B q can be computed as 



s 
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Proof. Set 



oo 



Fj(z)= I r z dN j (t) = Y, T -^z, 3 = 1,...,P, 

k=l I 



and 

F(z) = I t~ z dN(t). 



It is easy to verify that F{z) and the Fj(z) are analytic on the half-plane z > a and they 
are connected by the formula 

F(z) = f[F j (z). 
i=i 

Moreover, Fj{z) is analytic on 5?e z > r if j {ji, . . . ,j q }, whereas 

q z — a 

extend analytically to the same half-plane. Furthermore, these functions are of at most 
polynomial growth on any strip r < a < Sfte z < b. Thus, F has also at most polynomial 
growth on such strips and it is meromorphic on z > r, having a single pole at z = a of 
order v. The hypotheses from Aramaki's Tauberian theorem are therefore satisfied, and the 
result follows at once from it. □ 

2.3. Lower order terms in some special cases. When the {A^}^ =1 arise as = (cj(k— 
1) + bj)Pi , where Cj, bj, (3j are positive constants, and one assumes /3 p > /3 p _i > • • • > Pi > 0, 
it is possible to improve the asymptotic formula (|2.9p by giving lower order terms in the 
asymptotic expansion. Given c = (cx, C2, . . . , Cp), b = (pi, 62, . . . , b p ), (3 = (/3i, /?2, • • • , f3 p ) £ 
I?., we are interested in this subsection in the asymptotic behavior of the counting function 

(2.14) lA(A) = #[(h,...,k p )eW: (cifci + fti)*^ + b 2 f* . . . (c p k p + b p f* < a} . 

For the constants in our expansions, we shall need the Hurwitz zeta function [lj p. 251]. It 
is defined for fixed a > as 



00 ^ 

(2.15) C ( z;a ) = £___ ^ >L 

It is well-known that (|2.15p admits meromorphic continuation to the whole complex plane, 
with a simple pole at z = 1 with residue 1 (cf. [U p. 254] or [HI p. 348]). In particular, 
when a = 1 we recover Q(z) = ((z, 1), the Riemann zeta function. Using the Euler-Maclaurin 
summation formula [U \TT\ [21] , one easily deduces the following asymptotic formula 

(2.16) V - = ^ + ((s;a) + 0(A~ S ), A ->• 00, when < s and s ^ 1. 

*-r* (k + a) s 1 — s 

0<k<\ y ' 



WEYL ASYMPTOTICS FOR TENSOR PRODUCTS OF OPERATORS 9 

Observe that Proposition 12.11 immediately yields the dominant term in the asymptotic 
expansion of (|2.14p . 

We begin with the analysis of the case p = 2. The proof of the following lemma is inspired 
by the classical Dirichlet hyperbola method [U p. 57]. 

Lemma 2.6. Let = (pi,/3 2 ) be such that < Pi < /3 2 . Then, 

(2.18) Dl{\) = aN^M/c 2 ) xk ((^M A l O(A^). 



Proof. Since D^(X) = D^J^X/^c^c^ 2 )), where e = (1,1) and d = (61/01,62/02), we may 
assume that ci = c 2 = 1. For ease of writing, we set = D 13 '-. We have that 

4(A) = e 

(fci+bi)/ 3 l(fc2+6 2 ) /3 2<A 

fci+6 1 <A 1 /(^i+/3 2 ) ^ ^ / fc2+b2 < A i/(/3 1+/ 3 2 ) ^ / 

= A^I 1)/3l/( g 2 (A 1 /(ft+&) - ftj) + A^I 2j/32//3l (A 1 /( /3l+/32 ) - 6 2 ) - + o(A7TO), 

where Ij, s {x) = J2k<x + ^i)~ s ■ The asymptotic formula ([2.16P gives 

C(A/A; 61) + ^ — + 0(ATOw; 
pi - P2 

2 

1 PoX^+^2 1 

= A" 2 c(pV&; 61) + - R — ■=- + o(A"i+^ ), 

P2 — Pi 

and similarly 



Pi — P2 

The relation (|2.18p follows on combining the three previous asymptotic formulas. □ 
In general, we have: 

Proposition 2.7. lei /3 = (pi, . . . ,/3 p ) G M p 6e suc/i i/iai /3 P > pp-i > • • • > Pi > 0. Then, 
the counting function \2.1J$ has asymptotics 

(2.19) = X/^'^^ + 0(A' 9 i + '- ,+/3 p ), A^oc, 

i=i 

w/Mtre Aj = A jM$ = cj 1 c-MtoCiPv/P,; bu/c u ). 
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Remark 2.8. In (I2.19p . some of the terms may be absorbed by the error term, only those j 
such that 

(p - < p x + ■ ■ ■ + p p 

occur in the sum. Of course, this always holds for j = 1, 2; thus, at least, we always have two 
leading terms in (|2.19p . 

Proof. The case p = 2 is Lemma 12.61 Assume the result is valid for p — 1, we proceed to 
show (|2,19p by induction. As in Lemma 12.61 we may suppose that c\ = c 2 = • • • = Cp = 1. 

For simplicity, we write = D^-, . Set a = ^j=2^i' d = (^2> b^, . . . , b p ) 6 and 

ff= (&,...,&) ER p +\ Write 

D|(A) = Ii(A) + / 2 (A) + 0(A 1 /( a+ ' 31 )), 



where 



fei+6i<A 1 /( a +"i) 



/ A \ Vft i 

In(X)= V - I -\^+iHDl(X a/{a + Pl) ) 



n, p =2 ( fc J + fe J )' j < aq/(q+,3i) 

and 

D%X) = # {(fe 2 , . . . , fcp) G PP" 1 : (k 2 + 6 2 )^ . . . (hp + bp)*** < A} 
v 1 

= J^ijA^ +0(A £ ^), A^oo, 
i=2 

with ylj = \\ 2<VVJ Lj C{Pv/Pj] b v ) for j = 2, . . . ,p. If we combine the latter with (|2.16p . we 
conclude that the asymptotic behavior of ii(A) is 

m-±*i £ (F^M+ o ^- I,,< ° +fl0 > 

p / fl.\(a+(3,-)/(,8j(a+ft))\ 



i=2 
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Observe that C := /3f 1 / °° t' 1 ' 1 ^ Dj^dt is absolutely convergent. We then have 



J 2 (A) = A 1 ^ 1 / r l /^dD%t) - X^+h Di(\ a ^ a+ ^) 



\l//3i />oo 

= CA 1 ^ 1 - ±— / t- 1 - 1 'hD*J(t)dt 

PI J\<*/{°'+Pi) d 
P I .«.\(a+ft')/(/3,(a+/3i)) 

= cx 1 ^ 1 - V J ^ J + cka^- 1 )/^!)). 

Thus, we have shown (|2.19p except for C = 17^=2 C(Pv/Pi', b v ). But this fact follows by 
comparison with (|2.17p . The proof is complete. □ 

Remark 2.9. In connection with Proposition 12.71 Estrada and Kanwal have given an inter- 
esting distributional treatment of the asymptotic expansions of type (|2.19p , which often leads 
to improvements in the error term when interpreted in the distributional sense (cf. Sec. 
5.3]). 

3. Counting functions for tensor products of pseudo-differential operators 

We now apply results of Section [2] to the spectral asymptotics of the tensor products of 
pseudo-differential operators, and their perturbations. We shall mainly refer to operators in 
the Euclidean setting. Parallel results for operators on compact manifold will be outlined 
at the end. For the sake of completeness, we begin with a short survey of the classes of M. 
Shubin, cf. US E3 E2] . 

3.1. Globally elliptic pseudo-differential operators. Write z = (x,£) G M 2n and < 
z >= (1 + \z\ 2 ) 1 / 2 = (1 + \x\ 2 + |£| 2 ) 1/2 . One defines the class of symbols T™(R n ), m G R, 
< p < 1, as the set of all functions a G C°°(M. 2n ) satisfying, for all 7, 

(3.1) \d]a(z)\ < C 7 < z > m -"H, ze R 2n , 

with constants independent of z. The corresponding pseudo-differential operator is defined 
by Weyl quantization as 

(3.2) Pu(x) = a w u(x) = — ^ J j(*-vK (^±1, A u(y)dyd^. 

Note that if the symbol a is a polynomial in the £ variables, i.e. P in (|3.2|) is a partial 
differential operator, then the estimates (13. ip force a(z) to be a polynomial in the x— variables 
as well, i.e. P is a partial differential operator with polynomial coefficients. 

Let us introduce the global Sobolev spaces H s (W l ), s G N, Hilbert spaces with the norm 

(3.3) ||«||«= \\x a D^u\\ < 00. 

\a\+\/3\<s 
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By interpolation and duality the definition extends to s G R, and we have P| s i7 s (R n ) = 
<S(R n ), Ua-H" 8 ^") = S'(R n ). The immersion i : s H s -> H l is compact for s > t. If a G 
r™(R n ), then a™ : # S (R") -»• # fi - m (R n ) continuously for every s G R, hence a w : 5(R n ) -)■ 
5(R rt ), <S'(R n ) — > «S'(R n ). In the following we shall assume that for large \z\, 

(3.4) a{z) = a m (z) + a m _ p (z), 

where a m {tz) = t m a m (z),t > 0. We then say that a is globally elliptic if 

(3.5) a m {z) ^ for z ^ 0. 

Operators with globally elliptic symbol possess parametrix. Namely, there exists b G r~ m (R n ) 
such that a w b w =I + Ri and b w a w =I + R 2 , where R U R 2 : S'(M. n ) -> S(R n ). It follows that 
a w : i7 s (M n ) — >■ i/ s ~ m (IR n ) is a Fredholm operator and then eigenfunctions, i.e. solutions 
of a w u = 0, do not depend on s G R and belong to <S(R n ). Passing now to spectral theory, 
we assume that a G T 7 p n (W l ), m > 0, is real- valued and globally elliptic with a m (z) > 0, 
for z ^ 0. Then P = o^u : H m (R n ) ^ L 2 (R n ) is self-adjoint. The resolvent is compact 
and the spectrum is given by a sequence of real eigenvalues A& — > oo with finite multiplicity; 
the eigenfunctions belong to <S(R) and form an orthonormal basis. The spectral counting 
function iVp(A) = #{k : Aj. < A} behaves as 

(3.6) N P (X) = A\ 2n ' m + 0{\ a ), A^oo, 
for some a < 2n/m, with 



( 3 - 7 ) A = 7^ / dz - 

A sharp form of the remainder in (|3,6p can be obtained when a G r m (R n ) = r™(R n ) admits 
an asymptotic expansion in homogeneous terms a ~ X^fceN a m-2k- Then, with ^4 as before, 

(3.8) JV P (A) = AX 2n/m + 0(A 2(n ' 1)/m ), 

see, for example, Helffer [171 p. 175]. In the sequel, we shall assume that P is strictly positive, 
so that < Ai < A2 < • • • • For P as before, we may define the complex powers P z ,z G C. 
They are trace class operators if 5ftez < — 2n/m, and, by analytic continuation, we define the 
zeta function associated to P as 

00 

(3.9) Cp ( z ) = Tr(p-*)=J2*k Z - 

k=l 

3.2. Spectral asymptotics for tensor products. To give a precise functional frame to 
the results in the sequel, we shall introduce first the tensorized global Sobolev spaces. Write 
now Xj,yj G M. nj , Zj = (xj,yj) G R 2n ^, j = l,...,p, n = m + ... + n p , x = (xi,...,x p ), 
y = (yi,...,y p ) G R p , z = (zi,...,z p ) = (xi,yi, ...,x p ,y p ). For s = (s 1 ,...,s p ) G R p , we define 
the tensor product of Hilbert spaces 

p 

(3.10) H*(R n ) = (g)H s i(R n i). 

3=1 
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When the components of s are non-negative integers, from (|3.3p we recapture as norm 

(3.11) \\u\\g= Y, \\x ai ...x^D^...D^u\\. 

\aj\+\Pj\<*j 

We have f] g H g (R n ) = S(R n ) and \J s H g (R n ) = S'(R n ). The immersion i : H g {R n ) -> H\R n ) 
is compact if s > t, i.e. Sj > tj for j = 1, ...,p. 

As announced at the Introduction, we consider now Pj = aj in R n i , j = l,...,p, with 
real- valued symbol aj € r" 1 ^ (R n ^'), > 0,, and a mj {z) > for z ^ in (|3.5)) : we further 
define 

(3.12) P = P 1 ®...®P p , 

as operator P : H s (R n ) — > P s_m (IR n ), m = (mi, ...,m p ), for every s£ M p . In particular, we 
have P : Pf™(IR n ) -> L 2 (M n ) and P : <S(M n ) S(R n ), S'(R n ) -> S'(R n ). Moreover, P is 
self-adjoint and strictly positive, if the factors Pj are assumed to be strictly positive. 

If we denote by {X^}^ =1 the eigenvalues of Pj, according to the Introduction, the eigenval- 
ues of P are of the form A^ ...A^ and the eigenfunctions are tensor products of the respective 

eigenfunctions, hence they belong to 5(R n ). 

It is worth observing that P can be written in the pseudo-differential form (|3,2p with symbol 
a(z) = ai(zi)...a p (zp). However, the estimate (13. ip fails in general, and the considerations of 
Subsection 3.1 do not apply in this context. For the case p = 2, we address to [5], cf. [51 [26] . 
where a calculus was achieved in terms of vector-valued symbols. Here, to find an asymptotic 
expansion for iVp(A), we shall use (II. 2p in combination with the analysis of Section [2j In 
fact, from (13.6P and (13.7p . we have 

(3.13) Np j ~ AjX^I^ with Aj = — J— f dZj. 

Writing £p. for the zeta function of Pj, we inmediately obtain from Proposition l2.lt 

Theorem 3.1. Let P = P\ <g> ■ ■ ■ ® P p be as above and let a = m&Xj{2nj/mj}. Let fur- 
ther ji, ... ,j u be the indices such that a = 2nj q /m,j q , q = l,...,u. Then, P has spectral 
asymptotics 

(3.14) N P (X)= £ II ^t SX l)^ A ^°°' 

K-^ K<2 rip — 

where Aj is given by (|3.13p . 

We remark that the case p = 2, j/ = 1 or y = 2, of Theorem 13.11 also follows from the 
results of [6], see also [5]). 

As far as the reminder in (|3.14p concerns, from (|3.6p and Proposition 12.51 we obtain 

v-l 

(3.15) N P (\) = X a YC q log q X + 0(X r '), 

q=0 



(3-17) //;,,.,,. |< C ■ -''J' "T • I'r j = h-,P, 
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for some rj < a. The coefficient C v -\ is given by (|3. 14[) and the other constants C q , q = 
0, v — 2, are determined by (|2.12p . (|2.f3p . and the values of the derivatives or poles of the 
zeta functions Cpj( z ) &t z = a, j = I, . . . ,p. 

Willing sharp values of rj in the remainder, we further assume that a,- € r mj '(R nj ') with 
a,j ~ X]feeN a mj-2fc and we use (|3,8p . Proposition 12.41 yields. 

Theorem 3.2. Let P = P\X ... ® P p be as above. Assume that there is an index I € {1, ...,p} 
such that 1n\jmi > (3 = max,-^j{2rij /rrij} . Then 

(\ 2n i 

for any rj with r\ > max{/3,2(n; — l)/m/}. 

The following example shows that the exponent rj = /3 is sharp in (|3,16p . 

Example 3.3 (Tensorized Hermite operators). For tensor products of Hermite operators it 
is possible to detect lower order terms in the asymptotic expansion (|3.16p . Namely, let us 
fix = f3 p ) with j3\ < ... < f3 p ,c= (ci, c p ), b = (b\,...,b p ), p-tuples of positive real 

numbers, cf. Subsection 2.3, and consider 

|(-^+x|)-| 

so that for Cj = 1, bj = 1, we recapture Hj in (II. 3p of the Introduction. The eigenvalues of 

Hj lCj> bj, as one dimensional operator, are Ajj? = Cj(k — 1) + bj, k = 1,2,.... We then define 
the tensorized Hermite operator 

By Proposition 12.71 we have for the corresponding counting function 

P 

(3.19) N(X) = D^(\) = AjX 1 ^' + 0(X^+' 1+ ^) 

3=1 

with Aj as in Proposition ^. 71 In particular, for p = 2, Cj = 1, bj = 1, j = 1, 2, we obtain (|1.8|) 
of the Introduction. 

3.3. Asymptotics for lower order perturbations. For simplicity, we shall assume that 
the factors Pj in P = Pi®...®P p are partial differential operators with polynomial coefficients: 

(3.20) Pj= £ c^.a^g, ^-er. 

I«il+|^j[<mj 

As before, we assume that Pj is elliptic, with principal symbol 

(3.21) P®(x,£)= c a],^T^ >0 for (x„^/ (0,0), 

\a.j\+\pj\=mj 

self-adjoint and strictly positive, j = 1, ...,p. We shall study 

(3.22) A = P + R, 
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where R is a partial differential operator with polynomial coefficients having lower order with 
respect to P, in the sense that, writing a = (ati, ...,a p ),/3 = (f5\, ■■■,/3 p ) G N n ,n = m + ...+n p , 

(3.23) R= Capx^DP. 

\<*}\+\Pj\<mj 
3=1, -,P 

Note that each term of the sum in (13.231) can be regarded as a tensor product: 

x d D^ = x^D^ ®...®Xp p D% p p , 

hence for every s G W, 

A=P+R: H 3 {R n ) ->• H*-™(R n ). 

We shall first construct a parametrix for A. In absence of symbolic calculus, we shall use in 
the proof a direct argument. 

Proposition 3.4. For every fixed integer M > 0, we can find B : H s (R n ) — > H s+m (R n ) for 
every s = (si, s p ) G W, such that BA = I + S',AB = I + S", where S', S" : iP(R n ) -> 

H°+M(R n ), with M = (M, M). 
Proof. Consider 

p- 1 = p- 1 ® ... ® p- 1 : P ? (R") -> P" ff7 "(R n ). 
We have P"M = P -1 (P + P) = I — 5 with 

5 = -P- X P : P*(R n ) -> P*+ r (R n ). 



Define then 



We have 



M-l 

P = 5^'P- 1 : P*(R n ) -> P* +A? (R n ) 
3=0 



A/-1 

PA = ^ 5j(7 - 5) = / - S M , 

3=0 

where S" = -5 M : P^(M n ) -)■ P ?+M (R n ). It is easy to check that P is also a right parametrix. 

□ 

Corollary 3.5. The solution u G S'(R n ) of Au = f e S{R n ) belongs to S(R n ) 

Proof. If u G <S'(M n ), then u G H s (R n ) for some s. Taking P as in Proposition 13.4^ we obtain 

BAu = (I + S')u = Bf, 

hence, u = Bf - S'u. We have Bf G S{R n ) and S"u G i/ m7 (K n ). Since M in Proposition 
13.41 can be fixed as large as we want, we conclude u G 5(R n ). □ 



Corollary 3.6. The operator A : P*(R n ) H g -™(R n ) is Fredholm, for every fixed s G M n . 

Proof. Let us apply Proposition 13.41 with M = 1. Since the inclusion P s+1 (R n ) — > H s (R n ) is 
compact, the Fredholm property is proved. □ 
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Let us assume now that the operator A in (|3.22[) is self-adjoint. It follows from the preceding 
arguments that the resolvent is compact and the eigenfunctions belong to <S(R n ). Assume 
further that A is strictly positive; we write < Ai < A2 < ••• for its eigenvalues and Na for 
its spectral counting function. We give below an asymptotic formula for A&. In the sequel we 
write X to mean that / = 0(g) and g = 0(f) are both valid. 

Theorem 3.7. Let A = P + R in (|3.22p be as above. We use for P the notation of Theorem 
I3.il namely we write a = maxj{2rij /rrij} and we assume that a = 2nj/mj for v indices. We 
then have 

(3.24) \ k x fc 1 /«(Iog jfe)-^" 1 )/", k -> 00 
and 

(3.25) N A {\) x X a log"- 1 A, A -> 00. 
Proof. We have 

pu|| 2 = \\AP~ l Pu\\ 2 < Ci||Pn|| 2 with C x = \\AP~ l \\ 2 c{L2) . 
On the other hand, using Proposition 13.41 we may write I = BA — S' and thus 

||Pu|| 2 = \\P(BA - S')u\\ 2 < 2\\PBAu\\ 2 + 2||P,S / n|| 2 | < C 2 (\\Au\\ 2 + \\u\\ 2 ) 

with 

C 2 = 2max{||PP|| 2 (L2) ,||P5 / || 2 (L2) }, 

where (^^(^(l 2 ) < 00 if M in Proposition 13.41 is chosen sufficiently large. We now rewrite 
the preceding estimates as 

(A 2 u,u) < (dP 2 u,u) and (P 2 n, u) < (C 2 (A 2 + I)u, u). 

Using the classical max-min formula for the eigenvalues of A 2 , P 2 and denoting here ^ the 
eigenvalues of P, we deduce 

A 2 < Ci^l and m 2 <C 2 (A 2 + 1). 
Hence A^ X As a final step in the proof, we apply the following lemma. 
Lemma 3.8. If the sequence < fi\ < ^2 < ■■■> ~ * 00 > admits counting function 

N(n) ~ rfi a log s fi, [i — > 00, 

with r, a > and s > 0, then 

(3.26) ii k ~ (~y /a k 1/a (logk)- s / a , k^oo. 

The proof of this lemma is a simple combination of Proposition 4.6.4, page 198, and Lemma 
5.2.9, page 219, from [25J and it is therefore omitted. Since for the counting function Np(fx) 
we have from Theorem 13.11 

W P ( M )~rM a (log/i) ,/ - 1 

for a constant r, we deduce from Lemma 13.81 for the eigenvalues \x k of P the asymptotics 
(|3.26|) with s = v — 1. Hence (I3.24h follows. The asymptotic formula (I3.25h can be easily 
deduced from (|3.24|) . we leave details to the reader. □ 
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The rough asymptotics f)3.25|) can hopefully be improved, as suggested by the result from 
[6], which gives Na(X) ~ iVp(A) in the case p = 2. Furthermore, we expect formula (|3.15p is 
invariant under lower order perturbations. On the contrary, the precise asymptotics (|3.19p 
for tensorized Hermite operators should be lost, after addition of lower order terms. 

3.4. Pseudo-differential operators on closed manifolds. We now look at pseudo-differential 
operators on closed manifolds. Let M\,Mi, ■ ■ ■ ,M p be closed manifolds with dim My = rij. 
We consider elliptic self-adjoint pseudo-differential operators Pj on Mj of order rrij and prin- 
cipal symbol a mj (xj,£j) > for (xj,£j) € T*Mj \ (Mj x {0}),j = I,—, p. We denote by 
dxjdt^j the natural volume form on the cotangent bundle T*Mj. Under these circumstances, 
Hormander's theorem [321 Chap. Ill] gives us the asymptotic behavior of each counting 
function Np^X) of the eigenvalues {A^}^^ of Pj. In fact, 

Np j (A) = 1 = AjX n ^ m i + 0(A^- 1 )/" 1 J), A -> oo, 

where 



A-<A 



(3-27) ^ = 7^7 / <M6> 3 = h 



As usual, Cpj denotes the zeta function of the operator Pj. 

Proposition 12 . 41 directly gives the spectral asymptotics of the operator P = P\®P2®- ■ -®P p 
on the closed manifold M = M\ x Mi x • • • x M p of dimension dim M = n = n\ + n<i + • • • + n p , 
whenever one of the counting functions Np t dominates all the others. 

Theorem 3.9. Let Pj be elliptic self-adjoint strictly positive pseudo- differential operator as 
above, j = 1, . . . ,p,. Suppose that there is I G {1,2, ... ,p} such that ni/mi > nj/rrij for all 
j 7^ I. Then, the spectral counting function Np of the operator P = Pi (g> Pi ® ■ ■ ■ ® P p has 
asymptotics 



(3.28) Np(X)= Yl l=U,nCp J (^M))A" i/m! +0(AH, A 

aS 1 »a< 2 »...aS^<a V J#» 



oo, 



where A\ is given by (|3.27p and r satisfies max{(n; — l)/mi, max,^ rij/mj} < r < n\jmi. 

For the special case of the tensor product of two elliptic operators with one counting 
function dominating the other one, the error term in (|3.2p improves that from [51 Thrm. 3.2 
] for bisingular operators. 

We leave to the reader statements and proofs for the counterparts of Theorem 13.11 (|3.15p 
and Theorem 13.71 in the setting of closed manifolds. 
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